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ABSTRACT 

X-ray observations of several active galactic nuclei show prominent iron K- 
shell fluorescence lines that are sculpted due to special and general relativistic 
effects. These observations are important because they probe the space-time 
geometry close to distant black holes. However, the intrinsic distribution of Fe 
line strengths in the cosmos has never been determined. This uncertainty has 
contributed to the controversy surrounding the relativistic interpretation of the 
emission feature. Now, by making use of the latest multi-wavelength data, we 
show theoretical predictions of the cosmic density of relativistic Fe lines as a 
function of their equivalent width and line flux. We are able to show unequiv- 
ocally that the most common relativistic iron lines in the universe will be pro- 
duced by neutral iron fluorescence in Seyfert galaxies and have equivalent widths 
< 100 eV. Thus, the handful of very intense lines that have been discovered are 
just the bright end of a distribution of line strengths. In addition to validating 
the current observations, the predicted distributions can be used for planning 
future surveys of relativistic Fe hues. Finally, the predicted sky density of equiv- 
alent widths indicate that the X-ray source in AGNs can not, on average, lie on 
the axis of the black hole. 

Subject headings: accretion, accretion disks — galaxies: active — galaxies: nuclei 
— galaxies: Seyfert — X-rays: galaxies 



1. Introduction 



The X-ray spectrum of many active galactic nuclei (AGNs) exhibits a strong emission 
feature at an energy of 6.4 keV that is due to the fluorescence of weakly ionized iron in a dense 
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and r elatively cold medium (e.g., iNandra et al.lll989l : iPounds et al.lllQQOl : iNandra fc Pounds 



19941 ). The strength of the emission feature requires that the region responsible for the 
fluorescence subtend a bout half the sky as seen from the illuminating X-ray source (e.g., 
George fc Fabianlll99ll ). As the variability prop erties of AGN place the X-ray so urce within 



10 Schwarzschild radii from the black hole (e.g.. iGrandi et al.lll992l: lUttley 120071 ) , it is likely 
that relativistic effects may alter the observed line shape ( Fabian et al. 19891 : Laor 1991) 
allowing it to be a powerful probe of space-time curvature. Several examples of relativistically 
broadened Fe Kg l i nes have beeri detected in the spectra of AGNs over the last decade (e.g.. 



Tanaka et al.lll995l : iFabian et al.ll2002bl : iRevnolds fc NowakI 120031 : iMillerl 120071 : iFabian et al. 



20091 ). Th ese detections are leading to n i easurements of black hole spins and accretion disk 



dynamics (IBrenneman &: Reynolds! l2006l : [Reynolds fc Fabian! !2008! ) . 

The strength of the Fe lin e emission depends on only four parameter s: (i) the shape 

of the illuminating continuum (IBallantyne et al.! l2002l : iBallantyne fc Ross! !2002! ) , which in 



AGNs is a power-law with p hoton index F; (ii) the abundance of iron (relative to the solar 
value) in the accretio n disk (|Ballantvne et al.!!2002! ). Ap^', (hi) the ionization parameter of 



the illuminated disk (jBallantyne et al.!!20o3 i. ^ = AirFx/nu, where Fx is the illuminating 



flux, and nn is the density of the illuminated slab; and (iv) the reflection fraction R, a 
measure of the relative strength of the reflection component in the observed spectrum (this 
is related to the covering factor of the accretion disk). Recently, multi- wavelength survey data 
has been able to measure correlations between each of the flrst three of these parameters 
t o A = Lhn]/I/F,HH, the Eddi ngton ratio of AGNs (IRisaliti et al.! !2009! : Ilnoue et al.! 12007! : 



Netzer fc Trakhtenbrot! 12003). Here, Lboi is the bolometric luminosity of an AGN, and 



LEdd = 1-3 X 10^^ (Mbr/Mq) erg s ^ is the Eddington luminosity of a black hole with 
mass Mbh and Mq is the mass of the Sun. In addition, various surveys have been able 



to measure t 


le black hole 


(Labita et al. 


2009), as well 


Ueda et al. 


20031). We can 



2009!), as well as the luminosity function of AGNs as a function of redshift (e.g., 
sl). We can then combine all this information to calculate the density and 



evolution of the Fe Ka line from accretion disks. In the following we concentrate solely on the 
Fe Ka line that arises fr om the inner accretion disk, and neglect the contribution from any 
narrow compo nent (e.g.,lNandra!l2006l'). or one that arises from reflection off a Compton-thick 
absorber (e.g. iMurphv &: Yaqoob! 12009! ). The following A-dominated cosmology is assumed 
JSpergel et al1l2003h : Hn = 70 km s"^ Mpc'S Qa = 0.7, and Qm = 0.3. 
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Calculations 



The 2-10 keV X-ray luminosity function measured by lUeda et al.l (120031 ) is used as the 
basis for the density and evolution of AGNs. We consider 261 X-ray luminosities between 
10^^'^ and 10*^^ erg s~\ and assume a bolometric correction of 50 to convert these to Lboi 
(IVasudevan fc Fabian 2007 ). The black hole mass distribution at z ^ 0.15 measured by 
Netzei] J2009h is converted into 28 mass bins (spanning Mbh = 1-7 x 10^ Mq to 9.1 x 10*^ Mq) 
and normalized to sum to unity. Applying this distribution to each of the 261 luminosities 
yields 261 Eddington ratio distributions. At z = 0.15 the minimum and maximum Eddington 
ratios are A = 1.4 x 10~^ (at a X-ray luminosity of 10^^'^ erg s~^) and 2 x 10^ (at a X- 
ray luminosity of 10^^ erg s^^), respectively. A constant density ioni zed reflection code is 
used to calculat e the reflection contin uum including the Fe Ka line ( iRoss fc Fabianlll993 



Ross et al.lll999l : IRoss fc Fabianll2005l ). In these models, the density of the accretion disk 



is taken to be nn = 10^^ cm and Fx is changed in or der to obtain the desired value 



of £. The observed relationships b etween ^, F, Ap^ and A ( iRisaliti et al.ll2009l : llnoue et al. 



20071 : iNetzer &: TrakhtenbrotI 120071 ) are used to determine the input parameters for the the 
reflection code over a wide range of A. The parameters are frozen at reasonable valued for 
very high or low values of A. The relationships are: 



0.58(logA + 1.0) + 1.99 
2.55 



-2.0 < logA < -0.035 
log A > -0.035 



A 



Fe 



']^Q0.71ogA+0.59yQ_g 



10.45 



0.01 < A < 2 
A > 2 



(2) 



and 




(3) 



log A < -1 
-1 < logA < 1.5 
log A > 1.5. 

The Fe Ka equivalent width (EW) and line flux are measured for each A and then averaged 
over the A-distribution to obtain mean values for that particular X-ray lu minosity. The 
accretion disk is not expected to be optically thick to X-rays when A < 0.01 (INarayan Sz Yi 



^The upper lim it of F = 2.55 was chosen for computational reasons, but, despite the rare example (e.g., 
Fabian et al. 2009[ l. the vast majority of AGNs have photon indices less than this value (e.g., Risaliti et al 
2n09l ). 
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1995 ). and in those cases the Fe Ka hne EW and hne flux are set to zercj^. These calculations 
are repeated for 40 different redshifts betwee n z = 0.05 and z = 2, with the black hole mass 
distribution moving to larger masses as z^'^ ( iLabita et al.ll2009l ). 



It is currently unknown where the X-ray source is located in relation to the accretion 
disk in AGNs, and so we consider two limiting cases for R. T he first is to assume R = 1, 
a commonly observed value in measurements of broad Fe lines (INandra et al.l 120071 ). for all 
calculations of the Fe Ka EW and line flux. This is equivalent to assuming the accretion 
disk subtends a solid angle of ~ 2n as viewed from the X-ray source. The illuminating 
power-law and reflection spectra are simply added together to form the R = 1 spectra. The 
second case assumes that R is inversely proportional to A, which is a prediction of models of 
relativistic light bending in AGNs w here the X-ray source is located on the axis of the central 
black hole (IMiniutti fc Fabian! |2004| ). The relationship between R and A is determined from 
the Suzaku measuremen ts of the Seyfert 1 galaxy MCG-6-3 0-15 that assurned a variable 
R jMiniutti et al.ll2007h . A black hole mass of 4 x 10^ Mr;. dMcHardv et al.l 120051). and a 
2-10 keV bolometric correction of 50 ( Vasudevan fc Fabian 20071 ) was used to calculate A 
for the two states of MCG-6-30-15. A straight-line fit yielded: 



R 




A < 0.09 

0.09 < A < 0.92 

A > 0.92. 



(4) 



In both cases, the Fe Ka equivalent widths and line fluxes are computed by integrating 
over the total spectrum from 6 keV to 7.1 keV. For the EW calculation, the continuum 
underneath the line is estimated by fitting a straight line to the spectrum between those 
two energies. To calculate the line fl ux an emitting area from the disk had to be assumed. 
Following the line emission model of iNandra et al.l (120071 ). we assumed that each line was 
emitted from 1.235-6 with a flat emissivity index, and from 6-400 with an emissivity 
index of —3. Here Vg = GM-qy{/(? is the gravitational radius of a black hole with mass Mbh, 
G is the gravitational constant and c is the speed of light. 



^Ob servations of Galactic black hole binaries seem to indicate that inner edge of the disk is a function of 



A (e.g., 



Gierlinski et al. 


2008) 


e.g.. Miller et al. 


2006|) 



20081) . even when A > 0.01, although other observations imply a more stable inner 



2006h . If such an effect occurs in AGN disks, then the number density of relativistic 



lines from the lower A sources will be reduced. However, there is currently no observational evidence of such 
an effect in Seyfert galaxies or quasars. 
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Results 



At any multiplying the observed X-ray luminosity function by the slope of the com- 
puted luminosity-EW relationship results in the Fe Ka EW analogy to the AGN luminosity 
function, which we call the Fe Ka EW function. Figure 1 plots this function at four different 
redshifts. There are strong peaks in the i? = 1 function at EWs of ~ 110 eV and ~ 300 eV 
at all redshifts that correspond to strong neutral and ionized Fe Ka lines, respectively. In 
contrast, the strongest peak in the light-bending model is at an EW of ~ 500 eV which 
corresponds to neutral line emission that is enhanced by strong reflection. The R= \ model 
indicates that Fe Ka lines with EWs < 110 eV will dominate the space density at all red- 
shifts, while the light-bending model predicts a significant space density of very strong Fe Ka 
lines with EWs up to ~ 500 eV. 

Each panel in Fig. 1 also gives the average Fe Ka EW at the specified redshift for the 
two different models. The average EW for the i? = 1 model is ~ 50 eV for z < 1 and then 
rises to ~ 70 eV at higher redshifts. In the extreme light-bending model, the average Fe Ka 
EW is ~ 200 eV over the entire redshift range. 

Figure 1 also shows a clear difference in the expected density of relativistic Fe Ka 
with EWs > 100 eV between the R = 1 and the light-bending models. Although there 
exists up to a factor of 10 scatter in the observed relationships con necting F, ^, Ap^ and A 
(IRisaliti et al.ll2009l : llnoue et al.ll2007l : iNetzer fc TrakhtenbrotI 120071 ). this scatter could in no 
way explain the three order of magnitude difference between the R = 1 and light-bending 
models. Moreover, the decrease in high EW Fe Ka lines in t he -R = 1 mod e l is dr iven by the 



observed X-ray luminosity function. A recent compilation by lNandra et al.l (120071 ) found that 
the majority of broad Fe Ka components in nearby Seyferts have EWs less than 100 eV, in 
agreement with the Fe Ka function for R = 1. At z = 0.05, the light-bending model predicts 
a space density of AGNs with Fe Ka EWs > 100 eV that is 5x larger than the R = 1 model 
(this ratio drops to ~ 1 at 2; = 1); for EWs > 200 eV the space-density is 235 x larger in the 
light-bending model ( this ratio drops to y 150 at z = 1). Observations indicate that such 
intense lines are rare (INandra et al.l 120071 ). despite the fact that they would be the easiest 
to detect. However, our calculations assumed on ly one disk-like refle ctor in an AGN, but 
this may be violated for AGNs in certain states (jFabian et al. 2002a) where multiple disk 
reflections can boost the observed EW. Such a sce nario seems to be r equired to explain the 
very large Fe Ka EW observed from IH 0707-495 flFabian et allboOQl ). 



The predicted sky density distribution of relativistic Fe Ka lines in the range z = 0.05- 
0.5 is shown in Figure 2. This plot dramatically illustrates that, in the R = 1 scenario, 
AGNs with Fe Ka EWs > 200 eV will be very rarely observed between z = 0.05 and 0.5. 
On the other hand, if the X-ray source in AGNs is located above the black hole, as in the 
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extreme light-bending model, then such large EWs would be common, with sky densities 
of ~ 500 deg~^. If this model was generally true, then a near flat distribution of EWs 
should be obse rved, while, in con trast, the majority of detected broad Fe Ka lines have 



EWs< 100 eV (iNandra et al.l 120071 ). in agreement with the R = 1 curve. The bottom line is 
that, unless extreme light-bending is increasing the value of R, relativistic Fe Ka lines with 
EWs > 300 eV can only be produce by ionized accretion disks with large Fe abundances. 
Such a situation is most common when an AGN has a large Eddington ratio, and such AGNs 
are rare at z < 1. 

One can use these results to plan future Fe Ka line surveys performed by the proposed 
International X-ray Observatory (IXO). Figure 3 plots the predicted sky density of Fe Ka 
lines greater than a certain line flux in AGNs between z = 0.05 and 2. The plot shows 
that both the -R = 1 scenario and the extreme light-bending model predict that a small 
number of very bright lines with fluxes > 10~^^ erg cm~^ s~^. For reference, the Fe Ka 
line in the Seyfert 1 galaxy MC G-6-30-15 is one of the b rightest known and has a line 



flux of ~ 2 X 10 erg cm ^ s ^ ( IVaughan fc Edelsonll200ll ). This plot predicts that there 



should only be ~ 8 AGNs over the entire sky with Fe Ka lines that have fluxes equal to or 
greater than MCG-6-30-15, in agreement with the observed scarcity of such lines. At lower 
fluxes, the predictions of the two models diverge with the difference approaching 2 orders of 
magnitude at fluxes of ~ 10~^^ erg cm~^ s~^. Such a difference will be easily distinguished 
by future surveys by IXO, but could be measured by a careful analysis of the XMM-Newton 
and Suzaku archives. 



4. Concluding Remarks 

This work has shown, for the first time, the cosmic density and EW distribution of 
relativistic Fe Ka lines. Sensitive observations of these line profiles are vital as they allow 
measurements of fundamental parameters such as the spin of the central black hole. Our 
results validate the previous and current results from XMM-Newton a nd Suzaku that hav e 



shown that the majority of relativistic Fe Ka lines have EWs< 100 eV (INandra et al.ll2007l ). 
Now that models of the intrinsic line EW and flux distributions are available, detailed plan- 
ning of future Fe Ka surveys by IXO can be performed. The sensitivity of IXO will push 
spin measurements beyond the local universe and out to moderate redshifts, revolutionizing 
our understanding of the cosmic evolution of black holes. 

Finally, we have shown that the extreme light-bending model predicts many more intense 
relativistic Fe Ka lines than the R = 1 model (Figs. 1 and 2), in disagreement with current 
observational constraints. This implies that, on average, the X-ray source in AGNs does not 
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lie on the axis of the black hole where light-bending would be so extreme. It seems more 
likely that the X-ray source in average AGNs lies above the accretion disk at a radius of 
several Vg from the black hole where light bending is less severe. Comparing data from future 
surveys of Fe K« lines with plots such as Figure 3 should allow a measurement of the average 
radial displacement of the X-ray source. Of course, this conclusion does not preclude the 
possibility that extreme light-bending occurs in a few rare sources. 

The author thanks J.M. Miller and C.S. Reynolds for comments and advice. This paper 
is dedicated to the memory of Don Ballantyne. 
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Fig. 1. — The rest-frame Fe Ka equivalent width function ai z — 0.05, 0.50, 1.0 and 1.5. 
The sohd hne plots the function assuming a constant R — 1, while the dashed hne follows 
the strong hght-bending prescription. The peaks in the curves indicate EWs that are most 
common at that redshift. For the R = 1 model the peaks correspond to neutral Fe Ka lines 
(at ~ 110 eV) and ionized reflection (at ~ 300 eV). In the light-bending model, the lines 
are always produced by neutral iron. The roll-over at low EWs in the light-bending model 
at z = 0.05 corresponds to ionized Fe Ka lines from luminous quasars (which have very low 
space density at that redshift). The average EWs for both models are also indicated in the 
plots. The curves have been smoothed for clarity. 
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Fig. 2. — The expected sky density of relativistic Fc Ka lines from AGNs in the range z = 
0.05-0.5 with EWs greater than a certain observed value. The solid line plots the distribution 
from the models with a constant R = 1, while the dashed line shows the predictions from 
the light-bending model where R is inversely proportional to A. The EWs in this plot are 
measured in the observed frame. 
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log ^Fe Ka \-^^^ ^^'^ ^~^] 

Fig. 3. — The expected sky density of relativistic Fe Ka lines from AGNs in the range 
z — 0.05-2 with hne fluxes greater than a certain value. The sohd hne plots the distribution 
from the models with a constant R — 1, while the dashed line shows the predictions from 
the light-bending model where it! is inversely proportional to A. 



